The Fundamental Vibration of Molecular Hydrogen 
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The fundamental ground tone vibration of H2, HD, and D2 is determined to an accuracy of 
2 X lO"* cm~^ from Doppler-free laser spectroscopy in the coUisionless environment of a molecular 
beam. This rotationless vibrational splitting is derived from the combination difference between 
electronic excitation from the X^Ti^ ,v = and v = 1 levels to a common EF^'S^ ,v = level. 
Agreement within la between the experimental result and a full ab initio calculation provides a 
stringent test of quantum electrodynamics in a chemically-bound system. 



Quantum electrodynamics (QED), the fully quantized 
and relativistic version of electromagnetism, solves the 
problem of infinities associated with charged point-like 
particles and includes the effects of spontaneous particle- 
antiparticle generation from the vacuum. QED is tested 
to extreme precision by comparing values for the electro- 
magnetic coupling constant a obtained from measure- 
ments of the (7- factor of the electron [1 1 and from inter- 
fcrometric atomic recoil measurements 2| . These experi- 
ments and the Lamb shift measurements in atomic hydro- 
gen [3|,|4| have made QED the most accurately tested the- 
ory in physics. Concerning molecules, significant progress 
has been made recently in theoretical [5| and experimen- 
tal [6|, lZI investigations of QED phenomena in the HD+ 
molecular ion, where multiple angular momenta (rota- 
tional, electronic and nuclear spins) play a role. Neutral 
hydrogen has also recently been targeted for QED-tests, 
via a measurement of the dissociation energy of the H2 
[8|, HD [9|, and D2 [lO| molecules, and the experimental 
determination of rotationally excited quantum levels in 
H2 



11|. 



The rotationless fundamental ground tone (i.e. the vi- 
brational energy splitting between the v" = 0, J" = 
and v' = 1,J' = quantum states) of the neutral hydro- 
gen molecule is an ideal test system for several reasons. 
The total electronic angular momentum is zero for the 
X-^E+ ground state and the total nuclear spin for the ro- 
tationless J = state of para-H2 is also zero resulting in 
a simple spectrum without hyperfine splitting. The hy- 
perfine splitting is extremely small in HD (down to the 
Hz level [1^) and D2 in the absence of an / • J inter- 
action for the J = ground state. The recent progress 
in theory allows for calculations involving relativistic and 
QED-effects up to order a'* 13|,llJ]. Energy contributions 
in the calculation cancel to a large degree for the funda- 
mental ground tone, leading to a significant reduction in 
the uncertainty, thereby allowing for accurate QED tests. 

The present study focuses on a precise laser spec- 
troscopic measurement of the rotationless fundamental 
quantum of vibration in H2, HD and D2. In the ab- 
sence of rotation a one-photon transition between the 
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FIG. 1. (Color online) A schematic layout of the experimental 
setup. The oscillator cavity is seeded by a cw Ti:Sa laser, the 
pulsed output of which makes multiple passes in an amplifier 
stage. The amplified output is frequency up-converted in two 
frequency doubling (SHG) stages leading to fourth harmonic 
generation of ~211 nm. The deep UV radiation is sent to 
the experiment, where molecules in the X^E^ «' = 1 state, 
populated by electrical discharge, are optically excited in a 
two-photon Doppler-free configuration . The cw-seed light 
is compared to a frequency comb while the frequency offset 
between pulsed and cw-seed light is measured via on-line chirp 
analysis to obtain an absolute frequency calibration. See text 
for further details. 



(v" = 0, J" = 0) and (u' = 1, J' = 0) quantum states 
is strictly forbidden by nonrelativistic quantum mechan- 
ics. It is the nuclear spin-rotation coupling which makes 
this transition possible for D2 and HD, but the oscil- 
lator strength is extremely weak, while it remains zero 
for H2. Therefore an experimental approach is adopted 
to measure this quantity via the combination difference 
of separate two-photon transitions involving a common 
electronically excited state {EF^'S'^, z; = 0, J = 0) in a 
setup as displayed in Fig.[T] After a previous two-photon 
study of the Q(0) line in the EF-X{Ofi) band of the H2, 
HD and D2 isotopes 15|;ll6l, now the Q(0) two-photon 



lines in the EF — X(0,1) band are subjected to a fre- 
quency metrology experiment. The excitation scheme is 
indicated in Fig. [5] (b) . 

The experiment was performed with a long-pulse (20 ns 
duration) injection-seeded oscillator-amplifier Titanium 
Sapphire (Ti:Sa) laser system, specifically designed to 
deliver narrow bandwidth output (as low as 17 MHz) 
at its fundamental infrared wavelength [l7[. Its pulsed 
output was frequency up-converted in two frequency dou- 
bling stages to produce the resonant wavelength for in- 
ducing the two-photon transitions (211 nm for H2, 209 
nm for HD, and 207 nm for D2). The majority of the 
measurements were performed with a few /kJ per pulse 
of deep UV radiation, while up to 350 /xJ per pulse was 
used to assess the AC Stark effect. The deep UV radia- 
tion is arranged in a Sagnac interferometric configuration 
of counter-propagating probe laser beams, which facil- 
itates cancellation of the first-order Doppler shift |18| . 
The H2 is probed in the collisionless environment of a 
molecular beam, avoiding the usual difficulties of colli- 
sional frequency shifts in the weak electric quadrupole 
rovibrational spectra of hydrogen |19l421j . To reduce AC 
Stark effects a separate 355 nm laser, delayed by 30 ns 
with respect to the spectroscopy beam, is used to ionize 
the molecules in a 2-1-1' resonance-enhanced multi-photon 
ionization scheme. The molecular ions traverse a time of 
flight (TOF) tube before being detected. Population of 
the X^T,^, V = 1 state in the molecular beam is achieved 
through a pulsed pinhole discharge source [22 1 . 

Frequency calibration was carried out via a beatnote 
measurement of the Ti:Sa seed-frequency against a fiber- 
based frequency comb laser, resulting in an uncertainty 
smaller than 100 kHz. The dominant source of uncer- 
tainty is associated with the frequency offset, or chirp, 
between the cw-seed and pulsed output of the Ti:Sa laser. 
This phenomenon was previously characterized in detail 
for the setup [ITJ. For the present measurements the 
chirp was measured on-line for each laser shot, with av- 
erage chirp offsets of —5.00 MHz with an uncertainty of 
0.25 MHz. This leads to an uncertainty of 2 MHz on the 
EF — X transition frequencies, when a factor of eight 
is included to take into account the harmonic conversion 
and two-photon excitation. The systematic uncertainty 
due to the AC Stark effect was experimentally investi- 
gated by performing intensity dependent measurements 
of the transition frequencies and extrapolating to zero in- 
tensity (see Fig. [2]c). A listing of the uncertainty budget 
is given in Table HI 

A recording of the Q(0) two-photon line for H2 is shown 
in Fig. [21 (a) . The outcome of the present metrology ex- 
periments on the EF — X(0,1) transitions in the three 
hydrogen isotopomers is given in the lowest section of 
Table in combining this with results for the EF — X(0,0) 
band [15|, ll6| then yields the experimental values for the 
ground tone frequencies of H2 , HD and D2 listed in Ta- 
ble ini The experimental uncertainties are at the level of 



TABLE I. Estimated systematic and statistical uncertain- 
ties for the frequency calibrations of the EF^'E'^{v' = 0) - 
X^E+(u" = 1) transitions in H2, HD and D2 (in MHz). In 
the lowest section the Q(0) two-photon transition frequencies 
are listed in cm~^. 



Contribution 


Species Uncertainty (MHz) 


(i) AC-Stark^ 


H2 


0.5 




HD 


0.4 




D2 


0.8 


(h) DC-Stark 




<0.1 


(iii) Chirp'' 




2 


(iv) Frequency Calibration 




0.1 


(v) Residual Doppler 


H2 


1.1 




HD 


0.7 




D2 


0.6 


(vi) Statistics 


H2 


i.s-^ 




HD 


1.6"= 




D2 


2.2'* 


Combined Uncertainty'^ 




H2 


2.8 




HD 


2.7 




D2 


3.2 


Transitions 




(cm-i) 


EF -X{0,1) Q(0) 


H2 


95003.62059(10) 




HD 


95669.18608(10) 




D2 


96467.83195(10) 


EF-X{0,0) Q(0)' 


H2 


99164.78691(15) 




HD 


99301.34662(23) 




D2 


99461.44908(13) 



63 measurements. 
69 measurements. 
44 measurements. 

,(v)- 



^ Chirp and AC Stark offsets are corrected for and not in- 
dicated in the table. 
^ la statistical error based on 1 
'^ la statistical error based on 1 
"* la statistical error based on ■ 
'^ Quadrature sum of errors (i); 
* From previous studies [isl. [I6 



2 X 10^"* cm~^. This may be compared (see also Fig. [3]) 
with the values extrapolated from direct infrared spec- 
troscopy, yielding reasonable agreement with jl9|. How- 
ever there is a large disagreement at the 6cr-level with 
respect to the results of Bragg et al. [20[ , which were con- 
sidered the most accurate measurements to date. Good 
agreement is found with the less accurate results from 
Raman spectroscopy [23[, while the extrapolations from 
infrared measurements of HD [2J| and D2 [25| compare 
favourably with the present determination. 

To compare with theory, calculations of molecular hy- 
drogen level energies E are performed in the framework 
of nonrelativistic quantum electrodynamics (NRQED) 
making use of an evaluation in orders of the electromag- 



netic coupling constant a 

E{a) = [£(") + a^ f (2) ^ ^3 ^(3) ^ ^4 ^(4) 



(1) 



The nonrelativistic energy E^^' is obtained by first pro- 
ducing a Born-Oppenheimer potential with 15 digit ac- 
curacy [26,], and then solving the Schrodinger equation 
and calculating adiabatic and nonadiabatic corrections 
perturbatively in powers of the electron-nucleus mass ra- 
tio J27| . The resulting nonrelativistic binding energies are 
accurate to a few parts in lO"** cm~^ [27| and are in excel- 
lent agreement with the direct nonadiabatic calculations 
(a variational approach) for rotationless molecular hy- 
drogenic levels, performed by Adamowicz and co-workers 
for the case of H9 l28l| . For HD and D2 the agreement is 
less perfect 29l |30| , which we attribute to the need for 
much larger basis sets used in the direct nonadiabatic cal- 
culations {i.e. without relying on a Born-Oppenheimer 
approach), in particular for HD where the breaking of 
u ~ g inversion symmetry must be covered in the basis 
set. Moreover for heavier masses the calculations con- 
verge more slowly [3l|. The cancellation of errors for 
close-lying levels, probing the same part of the potential, 
leads to a significant improvement of uncertainties for the 
rotationless levels to 1 x 10~^ (see Table [Ill- 
Leading relativistic corrections S'^^^ are calculated from 
the expectation value of the Pauli Hamiltonian [13|, [iJ] . 
The leading QED corrections £'^^\ well-known in the hy- 



(a) 





0.5 


m 


0.4 






c 




3 


0.3 


n 








m 






0.2 


m 




c 




ai 


01 


U) 




c 











-0.1 




(b) EF'i; 



X'l 




-60 -40 -20 20 40 60 

f-28481 36895 MHz 



0.3 0.6 1 

Power Density (arb. units) 



FIG. 2. (Color online) (a) Recording of the EF'^EJ - 
X^E^(0,1) Q(0) two photon transition in H2 (in black). Be- 
low (in blue) the residuals from a fit to a Gaussian profile (red 
line) are shown, (b) Partial level scheme of molecular hydro- 
gen showing the measurement principle to obtain the ground 
tone vibrational splitting as a combination difference. The 
transitions measured in the present study are indicated by 
solid arrows while those measured in a separate study \XM, llQl 
are shown as dashed arrows, (c) Assessment of the AC-Stark 
effect in a plot of the absolute frequency of the EF — X(0,1) 
Q(0) transition in H2 vs. the 211 nm power density in the 
interaction zone. The open circle signifies the result shown in 
(a). 



drogen and helium atoms [32|, can be evaluated as ex- 
pectation values of more complicated operators, such as 
the Bethe logarithm [32| to an accuracy of 10^^ cm~^. 
The quoted uncertainties in £^^' originate from neglect- 
ing nonadiabatic and relativistic recoil corrections. The 
main issue at present is the S^'^'' QED correction, pro- 
ducing the largest contribution to the uncertainty. For 
molecular hydrogen this term has not been calculated ex- 
plicitly due to the high complexity of NRQED operators. 
Nevertheless, on the basis of results for the hydrogen and 
helium atoms, numerically leading contributions can be 
represented by calculable Dirac (5-functions. The remain- 
ing contribution, dominated by the electron self-energy, 
is represented by a correction to the Born-Oppenheimer 
potential, thus contributing to all molecular levels. A 
conservative uncertainty of 50% is estimated for contribu- 
tions related to the energy difference between vibrational 
levels. The radial nuclear functions for v — and v = 1 
probe almost the same range of internuclear distance, 
leading to significant cancellation of the uncertainty in 
S^'^' . This results in final theoretical predictions at an 
accuracy 1 x 10~^ cm~^ for the full QED evaluation of 
the rotationless fundamental ground tones in H2 , D2 and 
HD. All contributions to the energies and their uncertain- 
ties are presented in Table |ll] alongside the experimental 
results. 

The present experiment finds excellent agreement be- 
tween the measured values of the fundamental ground 
tones of the H2, HD, and D2 molecules and a full ab ini- 
tio calculation including nonadiabatic, relativistic and 
quantum electrodynamical effects at a combined preci- 
sion level of 2 X 10^^ cm~^. Implicitly included are cal- 
culations of electron correlations, a phenomenon known 
to pose a major difficulty in quantum ab initio calcula- 
tions of molecular structure. The comparison is graph- 
ically represented in Fig. [3l The present comparison is 
an improvement over the previous test of QED on the 
H2 dissociation energy [8[ (limited by the calculations to 
1.2 X 10^3 cm-i=36 MHz [H]) and over a test of high 
rotational states in H2 [ll| (limited by the experimental 
values at 5 X 10~^ cm~^=150 MHz). The result con- 
stitutes an accurate test of QED in a chemically-bound 
system. 

Although QED is considered to be the best tested 
theory to date, a recent disagreement on the proton 
size between muonic hydrogen and atomic hydrogen ex- 
periments opens a perspective on additional interac- 
tions between leptons and hadrons beyond the realms 
of QED 33, ISJ]. In view of this development, the high 
level of agreement between the most accurate theory and 
experiment for the molecular hydrogen level energies may 
be interpreted to constrain effects of possible long-range 
hadron-hadron interactions. Theories invoking extra 



compactified dimensions [35| , or predicting new particles 
and new interactions, of which a long-range interaction 
between hadrons is a possibility, may be parametrized by 
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FIG. 3. (Color online) (a) Current and previously measured 
values for the fundamental ground tone vibration in H2 [l9l . 
[20, l23l- The dashed horizontal line represent the value (and 
the yellow shaded area represents ±lcr uncertainties) from the 
full ah initio calculations including QED effects, (b) and (c) 
represent measured values for the ground tone of HD and D2 
respectively, along with the full ab-initio calculations. These 
are compared with the results of Rich et al. [241 ] and McKellar 
and Oka [25(| for HD and D2 respectively. 



a Yukawa-type potential V{r) = N'^^e^^^^^\ where N 
is the nucleon number, /? is a coupling constant and A is 
the effective range of the interaction. In the approxima- 
tion A ::^ _R = 0.74 A the present data yield a constraint 
on long-range hadron-hadron interactions, quantified by 
/3 < 6 X 10"^ eV-A. In this sense precision molecular 
spectroscopy opens an avenue to search for new physics. 
This research was supported by the Netherlands Foun- 
dation for Fundamental Research of Matter (FOM) 
through the program "Broken Mirrors & Drifting Con- 
stants" . We also acknowledge support by the NCN grants 
N-N204-015338 (J.K.) and 2012/04/A/ST2/00105 (K.P.) 
as well as by a computing grant from Poznan Supercom- 
puting and Networking Center, and by PL-Grid Infras- 
tructure. 



TABLE II. Summary of quantitative results with la uncer- 
tainties given in brackets (); all values in cm~^. Aq^*^ rep- 
resent the experimental values for the vibrational ground 
tones. The theory values are given for separate contribu- 
tions of Born-Oppenheimer energies Agi , adiabatic correc- 
tions Aof '"*" and nonadiabatic corrections Agi'^'^'^''^^ leading 
to a value for nonrelativistic energy /\Non-rei j^gjativistic ef- 
fects Agf ' and QED effects (leading Aq\ and higher order 
Aqj^ ) are added to form a final theoretical value for the 
ground tones: Aq^. Note that for each contribution, the can- 
cellation of uncertainties for the level energies of 11 = 1 and 
V — states is included. The bottom row lists the difference 
between experiment and theory for the fundamental ground 
tone vibration in the three hydrogen isotopomers. 



Experiment 



H2 



HD 



D2 



Kr 



4161.16632(18) 3632.16054(24) 2993.61713(17) 



Theory 



H2 



HD 



D2 



Non-relativistic contributions 
4163.40350 3633.71956 2994.44084 

-1.40284 -0.93259 -0.52150 

-0.83649 -0.62872 -0.30447 



aBO 

AAdiab 
01 

ANonadiab 
01 



A^ 



4161.16416(1) 3632.15826(1) 2993.61487(1) 



Relativistic and QED effects 
0.02341(3) 0.02093(2) 0.01771(2) 

-0.02129(2) -0.01863(2) -0.01539(2) 



Kf (a^) 



A««-° (a*) 


-0.00016(8) -0.00014(7) 


-0.00012(6) 


AS^ 


4161.16612(9) 3632.16041(8) 


2993.61708(7) 


Comparison 




H2 HD 


D2 


A Exp ATh 

^01 "^01 


0.00020(20) 0.00013(25) 


0.00005(18) 
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